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Abstract 

We consider binomial and inverse binomial sums at infinity and rewrite them 
in terms of a small set of constants, such as powers of vr or log(2). In order to 


perform these simplihcations, we view the series as specializations of generating 
series. For these generating series, we derive integral representations in terms of 
root-valued iterated integrals. Using substitutions, we express the iterated integrals 
as cyclotomic harmonic poly logarithms. Finally, by applying known relations among 
the cyclotomic harmonic polylogarithms, we derive expressions in terms of several 
constants. 

1 Introduction 

The goal of this article is to hnd and prove identities of the following form: 




-(l + 2bR0 81 x/3’ 

cxD ^ _2^ 
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Note that binomial and inverse binomial sums (see for example [5]) are of interest in 
physics: in particular, these sums have been studied in order to perform calculations of 
higher order corrections to scattering processes in particle physics [3, 19, 16, 13, 11, 12, 
15, 17, 21]. Central binomial sums were also considered in [18, 23, 9, 10], and there is a 
connection to Apery’s proof of the irrationality of C(3) (see [8]). In [22], a list of conjectures 
on series for powers of tt and other important constants is presented. In the frame of this 
article we are going to prove several of these conjectures involving binomial and inverse 
binomial sums. We summarize our approach with a concrete example. Consider the sum 


3 Si (fc) - 1 

h ^ 77 ) 


( 1 ) 


where 


S« (k) ^ 


2 = 1 


sign (a) ^ 

kM 


denotes the generalized harmonic numbers with a G Z \ {0}. As a hrst step we derive an 
integral representation for (1). Using tools from [5] this leads to 


3 Si (k) ^ 

h kHl^) 


'1 I rt I 


2 V ./q t Jq ^ s(4 — s) Jo 

L 1 C 


-drdsdt 


lo — t) Jo \/4 — s 


dsdt 


-1 


/o 1/^(4 — t) Jo 4 — s Jo \/4 — r 

r 1 


'I I rt 


t Jo \/s(4 - s) Jo y/r(4 - r) 


drdsdt + 1 

drdsdt. 


( 2 ) 


Next we want to rewrite these iterated integrals in terms of so called cyclotomic harmonic 
polylogarithms [2, 4, 6] which are extensions of the harmonic polylogarithms [20]. In 
order to dehne cyclotomic harmonic polylogarithms we introduce the following auxiliary 
function: For a G N and 5 G N, 5 < <p(a) (here (p denotes Euler’s totient function) we 
dehne : (0,1) i-G- M by 


if a = 6 = 0 
otherwise 

a{X) 

where $ 0 ( 2 ^) denotes the ath cyclotomic polynomial. 

Cyclotomic polylogarithms are now dehned recursively: Let m* = (a*, 6*) G bt < 
we dehne for x G (0,1) : 

H(x) = 1, 
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( ^(logx)^, if rrii = (0,0) for 1 ^ i ^ /c 

= S ^ (3) 

( lo faliy)'^m 2 ,...,mk{y)dy, Otherwise. 

The length k of the vector m = (mi, • • • ,mk) is called the weight of the cyclotomic 
harmonic polylogarithm Hm(x). If mi ^ (1,0) then the limit linij,^! Hmi,...,mu,(3:^) is hnite 
and we dehne 


hlmi,...,mu, (1) • Hmi(^) • 


H 


By nsing the snbstitntion x i+u+u'^ goiiig to show below how we can express 

the iterated integrals from eqnation (2) in terms of cyclotomic harmonic polylogarithms. 
Hence we get 



Finally, by using known relations (see [4, 6]) and new relations worked out in Section 4 
between cyclotomic harmonic polylogarithms at 1 we can derive 



(5) 


Summarizing, the proposed strategy, which has been implemented in the Mathematica 
package HarmonicSums^ [1], is 

1. Rewrite the sums in terms of nested integrals (see Section 2). 

2. Rewrite the integrals in terms of cyclotomic harmonic polylogarithms (see Section 3). 

3. Provide a sufficiently strong database to eliminate relations among these cyclotomic 
polylogarithms and hnd reduced expressions (see Section 4). 

Using this strategy the two main computational challanges are to transform the nested 
integrals into expressions in terms of cyclotomic harmonic polylogarithms, and to calcu¬ 
late the underlying relations such that the expressions in terms of cyclotomic harmonic 

^The package HarmonicSiims can be downloaded at 
http;//www.rise.jku.at/research/combinat/software/HarmonieSums. 
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polylogarithms collapse to simple constants. Here the computational effort grows ex¬ 
ponentially with the nested depth of the integrals. We are able express a large set of 
binomial and inverse binomial sums at infinity in terms of several constants using the 
proposed method. As a bonus, we can generate base identities that can be combined 
to new interesting results. Note that the main purpose of this article is to present this 
method which can be automated, hence not all identities presented in this paper are new 
identities. Many of the identities can already be found in [18, 23, 16, 13, 11, 12, 15, 17]. 
However we are also able to find new identies and prove several conjectures. In particular 
we are able to discover and prove some of the conjectures from [22]. 

The remainder of this article is organized as follows: In Section 2, we show in detail 
how we can derive integral representations for special binomial and inverse binomial sums. 
In Section 3, we show how these integral representations can be transformed to expressions 
in terms of cyclotomic harmonic polylogarithms. Section 4 deals with relations between 
the cyclotomic harmonic polylogarithms. Finally, in Section 5, we summarize the base 
identities that we found together with some nice combined results and list the conjectures 
from [22] that we could prove using our approach. 


2 Generating Functions and Infinite Nested Binomial Sums 


As a first step, we derive integral representations for the binomial sums. In order to 
accomplish this task, we view infinite sums as specializations of generating functions [5]. 
Namely, if we are given an integral representation of the generating function of a sequence, 
then we can obtain an integral representation for the infinite sum over that sequence if 
the limit x —?■ 1 can be carried out. This approach to infinite sums can be summarized 
by the following formula: 

OO OO 

i=l i=l 

In order to find integral representations of the generating functions, we rely on some 
results given in [5]. First we need the following well known properties. 


Lemma 1. Let IK be a field of characteristic 0 and /ef / : N —?■ K. Then the following 
identities hold in the ring IK[[x]] of formal power series 


n] = 


n=l 

OO 


E-/( 

n 

n=l 

OO n 

n=l 2=1 

OO Yi 

_ I iy _ I 0 


-| OO 

-Y,rf(n)dt, 


—E 


n=l 


n=l 


n=l 


x'" , 1 

n) - 

n ' X 


-f{n)dt. 


( 6 ) 

(7) 

( 8 ) 


n=l 


In addition, we make use of the following identities, which are useful for expressions 
involving binomial coefficients. Related formulae can also be found in the Appendix of 
[13], which do not explicitly express the results as iterated integrals. 
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Lemma 2. Let "K he a field of characteristic 0 and /et / : N —)■ K. Then the following 
identities hold in the ring IK[[a;]] of formal power series: 


E 


X 


n=l 


2n 


n 


E/w 




2 = 1 


Axl\- X >^0 


dt, 


u/i-t 


E;^E/w = E^/(") + 

n=l \ n ) i=l n=l \n) 


E 


X 




4 - X Jo - t) 

-dt. 


(2n +l)f") A'' ' A(4-x) Jo \/i(4-t) 

Let us illustrate the use of the formulae above by two simple examples. 
Example 3. Consider the generating function 


E 


X 


,=1 \ n J 


Applying (6) twice we obtain 


E 

71=1 


X 


1 ^ 


n 


t:) 


'0 


E 


n=l \ n J 


dt 


ri r 

L Jo 


S" 


E 


/ J /'2n'\ 

—f n( 

1=1 V 71 / 

Now, by virtue of (10) with f{n) := i, we obtain the result 


dsdt. 


E 


X 


1 


S "’h") 


t Jo ^/s{4: - s) Jo A(4^A 


drdsdt. 


Example 4. Consider the generating function 

x^'Si (n) 


E 


74=1 


n- 


A") ■ 


Applying (6) and then (10), we obtain 
x"Si(n) 


E 


1 ^ t^ 


71=1 


n‘ 


A") 


E 


dt A 


A 


E 


(9) 

dt, (10) 

(11) 


( 12 ) 


(13) 


(14) 


k Jo ofi^tJo A(4^fh"f:) 

Again, by applying (6) and then (10) with f{n) := 5n,i, we obtain 

(•tc 1 1 


dsdt. 


^ x”Si (n) 

2-^ 


71=1 


A”) 


t Jo s(4 — s) Jo \/4 — 


^drdsdt 
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+ 

+ 


r 1 x/i 

L ■\/t(4 — t) Jo \/4 — s 

/.XI /•* 1 , 




zdrdsdt + x I. 


/o ■\/t(4 — t) io 4: — s Jo \/4 — 

Note that combining Example 12 and 14 and setting x —)■ 1 we arrive at 


(15) 


1 - I Q (n) ^ 


a: 


E 


3 Si (k) - i 


1=1 \ n / n=l \ n / / k=l \ k J 

Since the limit a: —)■ 1 can be computed, we obtain (2) as our result. 


(16) 


3 Special Iterated Integrals to Cyclotomic Harmonic Polyloga¬ 
rithms 

In this section, we show how we can transform the iterated integrals that occur in the 
integral representation of Section 2 to expressions in terms of cyclotomic harmonic poly¬ 
logarithms dehned in (3). As an example, we consider the iterated integral 


lo Jo lo \/(4^^^(4 — s)v^4^^ 


drdsdt 


(17) 


which is part of (2). Substituting r —)■ in the innermost integral leads to 


lo Jo Ji 



— 2 — s+y/3\/ 4s — 
2(-l + s) 


\/3(l — u)^ 


y^(4 - t)t(4 - s){l + u + V?) 

Proceeding by the substitution s 

(1 -u) (- (\/3(l -u)2)) 


zdudsdt. 


(18) 


-2-f+v^V4f-f2 

r 2(-i+t) n 



0 Jl 


'i i/{A^^T)i {{1 + v) {1 + V + v^)) {1 + u + 


jdudvdt (19) 


and hnally by t ^ we arrive at 


0 rw nv 



(1 — u)(l — u)^ 


1 Jl Jl {I + w + w"^) {{1 + v) {1 + V + v"^)) {1 + u + 


rdudvdiu. 


( 20 ) 


Since the integration for cyclotomic harmonic polylogarithms always starts at 0, we rewrite 
this integral in the form 


■ dw 


1 — V 


dv I —y— "^J^du 


1 + w + w'^ Jq (1 u) (l-f u -1-u2) Jp (^i^u + u'^) 
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1 nv 


-3 

-3 

+3 


dw 



(1 — v){l — uY 



1 + W + w'^ Jo Jo ((1 + u) (1 + + t>2)) (1 + M + ^2) 

(l-uf 


Tdudv 


1 pW 


1 — V 


lo Jo (1 + w + ((1 + w) (1 + f 

'1 ( l _ y )( l _^)2 


dvdw 


■ du 



0 Jo (1 + W + M;2) ((1 + u) (1 + u + ti2)) (1 + M + ^2) 


0 (l + -u + -u2)^ 

- ^dudvdw (21) 


These integrals can be rewritten in terms of cyclotomic harmonic polylogarithms using 
partial fractions and the formula 




g{xy 


fib) 


/(a) 


* - 1 \9ib)^~^9'ib) g{a)^-^g'{a) 

fix)g"{x) f{x) 


Ja 9ix)^~^ V 9'ix)^ 9'ix) 

which can be derived by using integration by parts on 



dx 


( 22 ) 


Let us illustrate this on 

•*1 /•If 



la 9ix)^ ^9'ix) 


dx. 


lo Jo ((1 + n) (1 + w + n^)) (1 + M + 

First we look at the inner integral 


dudv. 


(23) 


{l-uY 


lo (1+M + M^)" 


rdu = 


lo (1+M + M^)" 

r 


jdu — 2 


u 


lo (l+M + M^)" 


rdti 


rdti 


Jo (i + u + u'^y 

Applying (22) to the first integral twice, together with partial fractioning yields 


(24) 


lo (l+M + M^)" 


rdw = 


{2v + 1 ) (^2 + y + 1 ) 
1 

{2v + 1 ) (^2 + y + 1 ) 


1 + 


1 + 


lo 3(2m + 1)2 3(^2+ M + 1) 


du 


Av 


3(2n + l) Jo 3(m2 + m + 1) 


du 


(-l + n)n ,2 

oH(3,o)('y). 


3 (1 + n + ^2) 3 

Applying the second integral and using the same strategy yields 


(25) 


(1 - v) 


(—1 + v)v 

lo (I + V {1 + V + v^)) \ 3(l+t> + n2) 


+ 2H(3^o)('?^)^ dv — --H( 2 ,o)( 1 ) + H(3^o)(1) 
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+ 2H(3,i)(1) + -H(2,o),(3,0)(1) - 2H(3,o).(3,0)(1) - 2H(3,i),(3,0)(1)- 

Proceeding in the same manner for the other two integrals in (24) and combining yields 
an expression in terms of cyclotomic harmonic polylogarithms for (23); 

4H(2,o),(3,0)(1) - 2H(3^o),(3,0)(1) - 4H(3_i)^(3^0)(1) 

-2H(2,o)(1) + 3H(3,o)(1) + 2H(3,i)(1) - 1. 

In a similar fashion we can rewrite all the iterated integrals in (21) in terms of cyclo¬ 
tomic harmonic polylogarithms. This leads to 

12H(3^0),(2,0),(3,0)(1) — 6H(3^o),(3,0),(3,0)(1) — 12H(3^o),(3,1),(3,0)(1) + (12H(3^o),(3,0)(1) 

+ 12H(3^0),(3,1)(1) + 12H(3^i)_(3^0)(1) - 12H(2,o),(3,0)(1) - 12H(3^o),(2,0)(1))H(3,o)(1) 

+ (12H(2,o)(1) — 12H(3^i)(1) — 6)H(3^o)(1)^ ~ 6H(3^o)(1)^ + 6H(3^o),(3,o)(1) + 1- (26) 

Note that the computational effort that has to be made grows exponentially with the 
nested depth of the integrals. In HarmonicSums sophisticated algorithms are used to 
speed up these computations. By applying these substitutions to all iterated integrals 
in (2) we can derive (4). 

Note that in our computations, integrands of the form \/4 — 3xy/x and \/2 — 
appear. In order to deal with these integrands, we used the substitutions x - 

and X 1 respectively. More precisely we utilize the following substitutions: 


x^ also 

-l+u^ 

-1—\^U+U^ 


y/(4 — x)x : 

X — )■ 

\/{A — 3x)x : 

X — >• 

a/(2 - x)x : 

X — )■ 

V^(l - x)x : 

X — )■ 

\/l — X : 

X — )■ 


1 + u + u'^ 
-1 + 

—1 — \/2u + 

{l-uf 

1 

1 

—Au 

(1 — uY 


4 Relations between Cyclotomic Harmonic Polylogarithms 

In this section, we examine relations beween cyclotomic harmonic polylogarithms at 1 
(compare [4, 6, 7]). A first set of relations originates from the shuffle algebra structure of 
cyclotomic harmonic polylogarithms i.e., the product of two cyclotomic harmonic poly¬ 
logarithms of the same argument can be expressed using the formula 

Hp(a;)Hg(a:) = ^ Hr(a;) 

r=p\AA q 
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(27) 
















in which piuq represents all merges of p and q in which the relative orders of the elements 
of p and q are preserved. As an example, we have: 

H(3,0),(2,0)(1)H(3,i),(3,o)(1) = H(3^o),(2,0),(3,1),(3,0)(1) + H(3,o),(3,1),(2,0),(3,0) (1) 

+H(3,0),(3,1),(3,0),(2,0)(1) + H(3,i),(3,o),(2,0),(3,0)(1) 
+2H(3_1)^(3^0),(3,0),(2,0)(1)- 

Secondly, there are the so-called dnality relations which are dne to argnment transforms 
of the form 1 — x—)-a;ora;—)■ |^, a, fe G M*. The transform x —)■ leads, for example, 

to the relation 


H(o, 0 ),(l, 0 )(l) ~ H(o,o),( 2 , 0 )( 1 ) + H(i^o),( 0 , 0 )(l) ~ H( 2 ,o),(l, 0 )(l) + H( 2 ,o),( 2 , 0 )(l) — 0. 

Other classes of relations for cyclotomic harmonic polylogarithm at 1 originate from the 
snm representation of these constants [6]. For a*, k E W, bi,n E N and q G Z* we dehne 

sign(ci)*i sign(cfc)*'' 




(^) = 




{aiii + {akik + bky^i^\' 


k is called the depth and w = Y^\^i\ci\ is called the weight of the cyclotomic harmonic 
snm S(ai,bi,ci),...,(afc,bi.,Cfc) ('^)- 

Fnrther relations can be discovered by ntilizing the connection between cyclotomic 
harmonic polylogarithms at 1 and cyclotomic harmonic snms at oo. This link can be 
established via the power series expansion of cyclotomic harmonic polylogarithms. In the 
following we sketch this translation mechanism, for details we refer the reader to [4, 6]. 
Let <ha(x) be a cyclotomic polynomial. Then we can write 


1 


a—1 oo 




Hence we get the power series expansion of depth one cyclotomic harmonic polylogarithms 
as follows. 


H(a,6)(a:) 


nx b rx °° 


y 


ai+q-\-h 


a—1 oo 


X 


q=0 
ai+q+b+l 


ai + q + b + 1 

ai+g—a+6+1 


^=0 2=0 
a—1 oo 

ai + q - a + b + 1' 

q=0 i=l 


X 


We can proceed recursively on the depth to obtain a power series expansion of a general 
cyclotomic harmonic polylogarithm. Let }lm{x) be a cyclotomic harmonic polylogarithm 
with 


H„(x) = 

j=i i=i 


^zi+Cj 

{zi + CjY^ ^ ) 
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for X G (0,1), w, gj G M, Cj G Z and some cyclotomic harmonic sums S„^. (i) . 
Provided that a = 0 or ak = z for some /c G M we have 


H( 0 , 0 ),m(^) — ^ ^ ^ ^ 


X 




-Sn, (*) , 


H 


(a,b),m 


{x) = 


^ ^ j,zi+q+b+Cj+l 

{zi + q + b + Cj + 1) 
j = l q=0 1=1 ^ j / 




where 


1 

*ha(a:) 


z—1 oo 

q=0 p=0 


Finally, for x ^ 1 these sums turn into cyclotomic harmonic sums at inhnity if c 7 ^ 1: 


E 




S.(^) 

{ai + by 


t S(a,b,c),n (^) • 


Example 5. Consider 

H(o,o),(3,2)(a:) 



CX) 

Ti = l 


(l + 3rOV’ 


Then, for a; —)■ 1 we get 

H(o,0),(3,2)(1) = gS(i,o, 2 ) ( 00 ) - S(34,2) ( 00 ) . 


In summary we can rewrite a cyclotomic harmonic polylogarithm at 1 in terms of 
cyclotomic harmonic sums at 00 using the power series expansion. From the sum repre¬ 
sentations, we can deduce additional classes of relations: cyclotomic harmonic sums form 
a quasi-shuffle algebra (see [4, 6 , 14]) and hence we get relations of the form 


S(5,3,3)(C)0) S(2,1,2),(3,2,1) ( 00 ) = 120 S(2,l,l),(3,2,l) (oo) -|- 19 S(2,l,2),(3,2,l) (oo) 

—45 S(2,i,2),(5,3,1)(oo) + 15 S(2,l,2),(5,3,2)(oo) — 5 S(2,l,2),(5,3,3) (oo) 

—300 S(5,3,i),(3,2,l)(c)o) — 100 S(5,3,2),(3,2,l)(oo) — 25 8(5,3,3),(3,2,1) (c)o) 

+ S( 2 , 1 , 2 ),(3, 2 , 1 ),(5,3,3) ( 00 ) S(2,l,2),(5,3,3),(3,2,l) (oo) 8 ( 5 , 3 , 3 ),( 2 , 1 , 2 ),(3, 2 , 1 ) (oo) . 


In addition to the quasi-shuffle relation there are two duplications: For q ^ 1 we have 


^ ^ 8(Q,^,b^,±c^),,,,,(Q,,^,6j,±ci) (2^.) 2 8(2am,bm,Cm)v,(2ai,6i,ci) ('^) 


where we sum on the left hand side over the 2 "* possible combinations. 
In addition, let di G { — 1,1}. Then we have 




dm ' ' ' diS(^am,bm.,dmCm),---,iai,bi,diCi) (2^.) 2 8(2a^,6^-am,Cm),.--,(2ai,f)i-ai,ci) ('^) 
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where we sum on the left hand side over the 2”^ possible combinations of the dt. 

Finally, there is the following multiple argument relation: For a, fc G N, 6 G No, c G Z*, 

k>2, 


k-l 



i=0 


and for a,, m G N, 6,, fc G Nq, c* G Z*, fc ^ 2, 



l,Cm-l),...,(ai,f)l,Cl) {k n) 

(fc j -i)sign(c™)^J-* 


{am{kj-i) + bi)\^\ 


As an example of a multiple argument relation, we state 


S( 2 . 1 , 2 ) (3n) = -S( 2 ,l, 2 ) (n) + S( 6 ,l, 2 ) (n) + S(6,5,2) (n) - ^ 


(6n + 5)2 225 ■ 


1 34 

, + 7 : 7 :; 


Note that if the lim^^oo Sm {n) exists then it makes no difference whether we consider 
sums at cxo or /c • cx), /c G N, and hence these relations remain valid in the limit. For 
example we get 


S(2,l,2) ( 00 ) — gS(2,l,2) ( 00 ) + S(6,l,2) ( 00 ) + S(6,5,2) (co) + 


Summarizing we get the following classes of relations [6]: 

• stuffle relations (quasi shuffle algebra of cyclotomic sums); 

• two duplication relations (if the sum is finite, it makes no difference whether the 
argument is cxo or 2 • cxo); 

• multiple argument relations (if the sum is finite, it makes no difference whether the 
argument is cxo or /c • cxo); 

• shuffle relations (shuffle algebra of cyclotomic polylogarithms); 

• duality relations of cyclotomic polylogarithms. 

As a first example, we consider cyclotomy 4, i.e., the cyclotomic polylogarithms with 


letters in 


{(0,0), (1,0), (2,0), (4,0), (4,1)} 


or the cyclotomic sums with letters in 
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w 

1 

2 

3 

4 

5 

^ sums 

4 

20 

100 

500 

2500 

# logs 

5 

25 

125 

625 

3125 

basis constants 

3 

1 

2 

3 

6 


Table 1: Number of cyclotomic harmonic polylogarithms and cyclotomic harmonic sums 
together with the number of basis constants of cyclotomy 4 at different weights. 

up to weight 5. Using the relations from above we can express all these constants using 
the following basis constants: 

H{i,o)(l), H( 2 ,o)( 1 ), H(4^o)(1); 

H(o,o),(4,i)(1); 

H(l,o),(0,0),(0,0)(l), H(o,o),(4,1),(4,0)(1); 

H(2,0),(1,0),(1,0),(1,0)(1); H(4^o),( 0,0),(0,0),(0,0) (1)) H(o,o),(4,l),(4,l),(4,0) (1)i 

H(i,0),{0,0),(0,0),{0,0),(0,0)(1),H(2,0),(1,0),(1,0),(1,0),(1,0)(1),H(4,1),(1,0),(1,0),(4,0),(1,0)(1), 

H(4,1),(1,0),(4,0),(1,0),(1,0)(1),H(4^1)^(4^0),(1,0),(1,0),(1,0)(1),H(4^1)^(4^1)^(1^0),(1,0),(1.0)(1)- 

As a second example we consider cyclotomy 6, i.e., the cyclotomic polylogarithms with 
letters in 

{(0,0), (1,0), (2,0), (3,0), (3,1), (6,0), (6,1)} 

or the cyclotomic sums with letters in 

[ (± 1 )^ (± 1 )^ (± 1 )^ 1 

1 i ’ 3i + 1 ’ 3i + 2 j 

up to weight 5. Using the relations from above we can express all these constants using 
the following basis constants (note that we only list the basis constants up to weight 4): 

H(i,o)(l), H( 2 ,o)( 1 ), H(34)(1), H(6,i)(1); 

H(o,0),(6,1)(1), H(6,1),(3,1)(1); 

H(o,0),(0,0),(1,0)(1)) H(o,o),(3,1),(6,1)(1)) H(o,o),(6,1),(3,1)(1)) H(o,o),(6,1),(6,0)(1)) H(o,o),( 6,1),(6,1)(1)) 
H(6,1),(6,1),(2,0)(1), H(e^i)^(6,i),(3,0)(l); 

H(o,0),(0,0),(0,0),(6,1)(1), H(o,0),(0,0),(6,1).(3,1)(1), H(o,0),(0,0),(6,1),(6,0)(1), H(o,0),(0,0),(6,1),(6,1)(1), 
H(0,0),(6,1),(3,1),(3,1)(1),H(0,0),(6,1),(3,1),(6,0)(1),H(0,0),(6,1),(3,1),(6,1)(1),H(0,0),(6,1),(6,0),(3,1)(1), 

H(o,0),(6,1),(6,0),(6,1)(1), H(o,o),(6,1),(6,1),(2,0)(1), H(o,o),( 6,1),(6,1),(3,0) (1), H(o,o),(6,l),(6,l),(3,l) (1), 
H(o,0),(6,1),(6,1),(6,0)(1), H(o, 0),(6,1),(6,1),(6,1)(1), H(2,0),(1,0),(1,0),(1,0) (1), H(6,l),(6,l),(3,l),(3,l) (1), 
H(6,l),(6,l),(6,l),(2,0)(l),H(6,i),(6,l),(6,l),(3,l)(l)- 

In the following we give some example relations for the cyclotomic harmonic polylog¬ 
arithms appearing in (4): 

H(3,0)(1) = H(6,i)(1), 
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w 

1 

2 

3 

4 

5 

# sums 

6 

42 

294 

2058 

14406 

# logs 

7 

49 

343 

2401 

16807 

k basis constants 

4 

2 

7 

18 

52 


Table 2: Number of cyclotomic harmonic polylogarithms and cyclotomic harmonic sums 
together with the number of basis constants of cyclotomy 6 at different weights. 


H(2,0),(3,0)(1) 

H(3,0),(1,0)(1) 

H(3,0),(2,0)(1) 

H(3,0),(3,0)(1) 

H(3,0),(3,1)(1) 

H(3,1),(3,0)(1) 

H(3,0),(1,0),(3,0)(1) 

H{3,0),(2,0),(3,0)(1) 

H(3,0),(3,0),(1,0)(1) 

H(3,0),(3,0),(3,0)(1) 

H(3,0),(3,1),(3,0)(1) 

H(3,1),(3,0),(3,0)(1) 


H(2,0)(1)H(6,1)(1) - - -H(o,o),(6,1)(1), 

- in,8.1,(1)2 - jH,8,1,(1), 

+ -H(o,o),(6,1)(1), 

2^(6, 

H(3,1)(1)H(6,i)(1) - -H(o,o),(6,1)(1), 

2 H( 0 , 0 ),( 6 ,i)( 1 ), 

12 2 

“2^(6,- 2H(6,i)(1)H(0,0),(6,i)(1) - ^H(0,0),(0,0),(1,0)(1), 

7 7 13 

-^H(6,i)(1)^ - -H(6_i)(1)H(0,0),(6,1)(1) - Y^H(o,o),(0,0),(1,0)(1), 

2H(3,i)(1)H(6,i)(1)^ - -H(6,i)(1)H(0,0),(6,i)(1) + ^H(0,0),(0,0),(1,0)(1), 

gH( 6 ,i)(l)^, 

4 5 14 

- 2 H( 6 ,i)( 1 )^ - -H(6,i)(1)H(o,o),(6,i)( 1) - ^H(o,o),(o,o),(i,o)(l), 

2 7 

2 H( 6 ,i)( 1 )^ + H(6,i)(1)H(o,o),(6,i)( 1) + ^H(o,o),(o,o),(i,o)(l)- 


Plugging these relations into (4) we find 


3 Si (/c) I 

k H7) 


~H(0,0),(0,0),(1,0)(1) — S(l,0,3) (oo) — C 3 . 


5 Results 

In this section, we provide a number of base identities we discovered using the techniques 
outlined in Sections 2, 3, and 4. These results were obtained and proved using our 
implementation in the Mathematica package HarmonicSums [1]. Note that this section 
should illustrate what can be achieved using our implementation. Many of these base 
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identities can already be fonnd in [18, 23, 16, 13, 11, 12, 15, 17]. We define here a nnniber 
of constants that appear. Note that these constants do not possess any fnrther relations 
indnced by the algebraic properties given in Section 4. 
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For weight 5 and 6 we just list a few nice base identities 
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In the following we state some identities that can be derived from these base identities. 
Combining (64) and (65) leads for instance to: 
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In addition, we were able to prove the following conjectures from [22]: 
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h *7?) 

= 5c3. 
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^ Si {2k) + 17Si {k) 

h ^??) 

= ttJ’C, 

(119) 

°° 2 ^ 

S *’??) 

(^2Si ( 2 fc) - 3Si {k) + 0 

= lC3. 

( 120 ) 

°° nk / 

6Si(2fc)-llSi + 

= 27rC, 

( 121 ) 

°° 2 ^ 

S *’??) 

(^2Si ( 2 fc) - 7Si {k) + 0 

TT^ 

= -ylog 2 . 

( 122 ) 

3 ^ 

S *’??) 

(^ 6 Si {2k) - 8 S 1 (fc) + ^ 

^fC3, 

(123) 

°° qfc / 

S fc??) ( 

6 S 1 {2k) - lOSi {k) + 0 

= 2\/3 7rJ’S:, 

(124) 


°° ok / 1 \ 

S *■(■;) (=■“’*-) 

ir^, 

= ylog3. 
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^ Cs + ^^vrCs, 
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Note that here L ;= 2c8 — ||g + 2 and K := 2ci — ^ + 2. 


(133) 

(134) 

(135) 
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